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Introduction We introduce a new tool to planetary ge- 
ology for quantifying the spatial arrangement of vent 
fields and volcanic provinces using non-parametric ker- 
nel density estimation. Unlike parametric methods where 
spatial density, and thus the spatial arrangement of vol- 
canic vents, is simplified to fit a standard statistical dis- 
tribution, non-parametric methods offer more objective 
and data-driven techniques to characterize volcanic vent 
fields. This method is applied to Syria Planum volcanic 
vent catalog data as well as catalog data for a vent field 
south of Pavonis Mons. The spatial densities are com- 
pared to terrestrial volcanic fields. 

Kernel Density Estimation Spatial density can be es- 
timated using a kernel function, which is simply a spa- 
tial density about a point, the center of a volcanic vent. 
These kernels are summed to calculate the spatial density 
of vents for any location in the region of interest. Kernel 
functions are dependent on a smoothing parameter, or 
bandwidth, which is analogous to standard deviation of a 
Gaussian distribution. For one kernel function, given by 
Wand and Jones [1], the bandwidth is given as a 2 x 2 
element matrix, H, used to estimate spatial density, A(s), 
at point s: 
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into volcanic fields [5,6]. Two previously cataloged 
shield fields include Syria Planum [6] and a field of 
volcanic vents south of Pavonis Mons [5]. In both of 
these projects, volcanic vents are defined as positive to- 
pographic features (10s to 100s m in relief) that exhibit 
radial flow features and craters at the apex of the feature. 
Vents are cataloged as 2-dimensional points to represent 
4-dimensional volcanic event that created the feature. 
The Syria Planum catalog contains N = 263 vents. The 
catalog of vents south of Pavonis Mons contains N = 96 
vents. 

Syria Planum and Pavonis Mons catalogs are each 
used to estimate H using the SAMSE bandwidth selec- 
tor. The square root of the bandwidth matrix for Syria 
Planum vents is 
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This matrix gives a E-W smoothing distance of 103.46 
km (twice the upper left element, 51.73) and a N-S 
smoothing distance of 95.9 km (twice the lower right el- 
ement, 47.95). The negative values in the upper right 
and lower left elements rotate the kernel ellipse counter- 
clockwise. The square root matrix for the vent field south 
of Pavonis Mons is 


\/H = 


24.43 0.31 

0.31 49.91 
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where |H| is the determinant of the bandwidth matrix, 
b = H -1 / 2 d and b T is the transform of b. d is a 1 x 2 
distance matrix, and N is the total number of volcanic 
vents. 

The bandwidth is optimized using a smoothed 
asymptotic mean integrated squared error (SAMSE) 
method, developed by Duong and Hazelton [2] and freely 
available in the statistical package, R. 

Tharsis Low-Shield Volcanic Fields Many new small 
shield volcanoes (10s km in diameter) have been identi- 
fied within the Tharsis Volcanic Province through post- 
Viking images and topography [3-6]. Using topographic 
data from the Mars Orbiter Laser Altimeter (MOLA) 
gridded dataset as well as images from the Thermal 
Emission Imaging System (THEMIS), the Context Im- 
ager (CTX), and the High Resolution Stereo Camera 
(HRSC), workers have begun cataloging volcanic vents 


giving a E-W smoothing distance of 48.86 km and a N-S 
smoothing distance of 99.82 km, with the elliptical ker- 
nel being rotated clockwise. 

These bandwidth matrices are used to map the spatial 
density of each vent field. These spatial density distribu- 
tions are mapped in Figures 1 and 2, which display the 
25 th , 50 th , 15 th , 95 th , and 99 th percentile of vent den- 
sity as contours. Note that these maps integrate to one. 
Spatial intensity (vents/km 2 ) can be estimated by multi- 
plying results of equation 1 by N. 

Discussion/Conclusion SAMSE bandwidth estimates 
for Syria Planum and Pavonis Mons vent catalogs yield 
smoothing distances within 10% of each other, with 
the semi-major axes of each elliptical kernel having a 
smoothing distance of ~ 100 km. The most striking fea- 
ture of these bandwidths is they are 1-2 orders of magni- 
tude greater than bandwidths for monogenetic volcanic 
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fields on Earth [7,8]. For example, vent density in the 
sparse Yucca Mountain, Nevada volcanic field results in 
a semi-major axis smoothing distance of 28.8km [7], 
Similarly, volcano density is greater among terrestrial 
arc poly genetic volcanoes. In the Tohoku arc, Honshu, 
Japan, we find a semi-major axis smoothing distance of 
~33 km. Thus, volcano spatial density in Syria Planum 
and Pavonis Mons catalogs is demonstrably lower than 
in terrestrial monogenetic volcanic fields and arcs. 

The resulting spatial density maps for Syria Planum 
and Pavonis Mons catalogs reveal single modes in vent 
density in each field (Figures 1 and 2). The areas of 
these volcanic fields can be estimated from the area en- 
closed by the 99 th percentile isoline for each spatial den- 
sity model. For the Syria Planum catalog this area is 
445,100 km 2 and for the catalog of vents south of Pavo- 
nis it is 160,700 km 2 . Again, in both cases these areas are 
greater than for terrestrial monogenetic volcanic fields or 
arcs. 

The SAMSE bandwidth estimates provide a means 
of objectively assessing anisotropy in vent distribution 
that may reveal tectonic processes. The Syria Planum 
kernel is roughly circular. The vent field south of Pavo- 
nis Mons has an elliptical kernel. Very low deviatoric 
stress is expected to produce a radially symmetric kernel. 
While both density distributions are unimodal in nature, 
perhaps implying that each was created by one uniform 
melt region, the elongation of the resulting Syria Planum 
density distribution and the elongated kernel of the vent 
field south of Pavonis Mons pose questions. Was there 
preferential dike alignment in the NNE-SSW direction 
south of Pavonis Mons? Was the underlying melt region 
of Syria Planum elongated or was the shape of the shield 
field created by near-surface rheology? Or do these ob- 
servations result from mapping biases or potential burial 
of sections of the Pavonis Mons field? 

Kernel density estimation is a new method to provide 
objective estimates of spatial density of vents and area of 
planetary volcanic fields, simplifying comparisons. We 
find Syria Planum and Pavonis Mons catalogs to indicate 
lower vent densities and larger areas than observed for 
terrestial volcanic fields, but are nevertheless character- 
ized by unimodal distributions. 
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Figure 1: Spatial density distribution of vents on Syria 
Planum. In both figures, contours represent the 25 , 
50 th , 15 th , 95 th , and 99 th percentiles of vent density; 
black circles represent volcanic vents; north is up; east- 
ings and northings are given in Transverse Mercator kilo- 
meters. 



Figure 2: Spatial density distribution of vents in a vent 
field south of Pavonis Mons. 
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